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We investigate magnetic properties of liglitly doped antiferromagnetic Mott insulators in the 
presence of non-magnetic impurities. Within the framework of the t-J model we calculate the doping 
dependence of the antiferromagnetic order parameter using self-consistent diagrammatic techniques. 
We show that in the presence of non-magnetic impurities the antiferromagnetic order is more robust 
against hole doping in comparison with the impurity- free host, implying that magnetic order can 
re-appear upon Zn doping into lightly hole-doped cuprates. We argue that this is primarily due to 
the loss of coherence and reduced mobility of the hole quasiparticles caused by impurity scattering. 
These results are consistent with experimental data on Zn-doped La2-i:Sra;Cu04. 



I. INTRODUCTION 

The undoped parent compounds of the high-Tc su- 
perconducting cuprates, like La2Cu04, are antiferromag- 
netic (AF) insulators. Upon introducing a small amount 
of holes, e.g., by substituting La with Sr in La2Cu04, 
the commensurate long-range order is rapidly suppressed 
and completely destroyed at few percent of hole doping 
(xc w 2% in La2_a;Sr^Cu04).^ For 2 % < a; < 5.5%, 
La2-2;Sr2;Cu04 enters a spin-glass phase with incommen- 
surate spin correlations, and finally becomes supercon- 
ducting for larger x, with only short-range AF correla- 
tions in the superconducting state. It is believed that a 
proper understanding of the magnetism is crucial for a 
theoretical description of high-Tc cuprates. 

On the theoretical side, the suppression of AF order 
upon hole doping is reasonably well understood. Nu- 
merous approaches employ effective single-band models 
with strong correlations to describe the copper oxide 
planes. The magnetic order, present in the half-filled 
Neel state, is scrambled by the motion of holes - this 
leads to heavy hole quasiparticles (so-called magnetic or 
spin polarons) and to a reduction of the magnetic or- 
der parameter, accompanied by a softening of the spin- 
wave spectrum. (A different approach describing the 
destruction of AF order, based on the idea that holes 
introduce strongly ferromagnetic frustrating bonds, has 
I been proposed in Refs. 6,7.) 

Localized impurities provide an interesting tool to in- 
' vestigate bulk properties of strongly correlated systems 
in general, and high-Tc cuprates in particular. Dop- 
ing of non-magnetic impurities (i.e. vacancies) into the 
parent antiferromagnet,^"^" e.g. in La2Cui_zZn204,^^ 
leads to a gradual suppression of magnetism, with long- 
range order surviving up to the percolation threshold, 
Zp « 40.5 %. One expects that combining hole dop- 
ing and non-magnetic impurities would lead to a de- 
struction of the AF order even faster than with hole 
doping alone. Surprisingly, magnetization measure- 
ments on La2-a;Sr2,Cui_zZnz04 by Hiicker et al}^ have 
shown that this is not the case. They found that in- 



troducing non-magnetic Zn impurities in lightly doped 
La2-2;Sr2;Cu04 (x < 4%) causes an initial increase of 
the Neel temperature, Tn, while only an addition of Zn 
beyond z « 10 % leads to a decrease of Tn. Concomi- 
tantly, the decrease of Tn with increasing hole doping x 
is much weaker in Zn-doped compounds. Recently, it has 
been found that such an impurity-induced re-appearance 
of the AF ordering is even more pronounced in samples 
with magnetic Ni impurities. ^'^ 

The purpose of the present paper is a theoretical de- 
scription of the combined effect of mobile holes and static 
vacancies in lightly doped Mott insulators. We argue be- 
low that vacancies primarily prevent the holes from de- 
stroying the long-range ordered magnetism; other (purely 
magnetic) effects are subdominant. Starting from a mag- 
netically ordered state, we consider the influence of mo- 
bile holes on the host magnetism in a situation where 
hole quasiparticles scatter off the static impurities. Tech- 
nically, we employ a, t-J model and the self-consistent 
Born approximation (SCBA) to account for the interac- 
tion between holes and spin waves;^'^ the effect of the 
impurities is treated via a self-consistent T-matrix ap- 
proach. We show that the destructive effect of holes 
on the background magnetism is reduced (i) when holes 
are less mobile or less dispersive and (ii) when holes are 
less coherent, both changes caused by the addition of 
non-magnetic impurities. This results in commensurate 
antiferromagnetism being more robust against hole dop- 
ing for impurity-doped host systems compared to clean 
ones. We obtain detailed magnetization curves as func- 
tion of hole and impurity doping, being in good agree- 
ment with the available data on La2_2:Sr2,Cui_2Znz04. 
(Incommensurate magnetism, as occurs for hole doping 
X > 4%, with much smaller ordering temperature, is not 
the subject of this paper.) 

Before describing our theoretical approach, we briefly 
summarize impurity effects in insulating quantum mag- 
nets and in cuprate materials. 
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A. Impurities in quantum magnets 

In a quantum spin system (i.e. a Mott insulator), which 
can be tuned between an antiferromagnetic and a para- 
magnetic ground state, the effect of static non-magnetic 
impurities (i.e. vacancies) is markedly different in both 
phases, the origin being the random Berry phases intro- 
duced by the doping. As mentioned above, introducing 
vacancies into the long-range ordered magnet weakens 
the magnetic order. At small impurity concentrations, 
this effect can be well captured using linear spin-wave 
theory.^ For a square lattice geometry, long-range mag- 
netism persists up to the percolation threshold^'* Zp - this 
has been nicely verified for La2Cui_2Zn204 using neu- 
tron scattering. The quantum phase transition at Zp 
displays an interesting interplay of classical percolation 
physics and quantum effects. 

In contrast, removing spins in the quantum paramag- 
net (which is a valence bond crystal) has been shown 
to create effective magnetic moments in the vicinity of 
the dopants. ^^"^^ Microscopically, the formation of the 
moments can be understood as breaking of host valence 
bonds. The impurity-induced moments are magnetically 
coupled via the (gapped) host spin excitations and - in 
the absence of geometric frustration - eventually order 
at very low temperatures. Thus, at zero temperature 
the paramagnetic phase of the clean host is replaced by 
a phase with weak magnetic order through doping with 
non- magnetic impurities. 

Besides vacancies, also magnetic impurities may be in- 
troduced into quantum magnets, either by replacing host 
ions by ions of different spin (e.g. Cu by Ni), or inserting 
additional (e.g. out-of-plane) magnetic ions. At small 
doping, the physics is similar to the one in the case of 
vacancies, i.e., it is dominated by the low-energy behav- 
ior of the impurity moments (random Berry phases), but 
percolation physics at larger doping is absent. 

B. Impurities in cuprates 

Impurity effects in cuprate superconductors are di- 
verse, most prominent being the rapid suppression of 
superconductivity by in-plane impurities like Zn or Ni 
substituting for Cu. For hole concentrations in the su- 
perconducting regime, it has also been shown that nom- 
inally non-magnetic impurities like Zn or Li induce mag- 
netic moments, which have a staggered spatial structure 
around the actual impurity site. NMR^^ shows that the 
impurity moments form at temperatures above 400 K, 
and the formation mechanism is believed to be similar to 
that in an insulating paramagnet, i.e., breaking of sin- 
glet bonds described above. However, impurity-induced 
magnetic order has not been detected, the reason likely 
being strong quantum effects (e.g. Kondo screening of 
the moments) due to the presence of mobile charge car- 
riers. (One exception is magnetic order in Co-doped 
YBa2Cu306+5^^ - due to the large spin of Co quantum 



effects arc weaker.) Further effects of Zn and similar 
impurities in cuprates with larger hole doping include 
magnetic in-gap states, ^^'^^ the broadening of the gapped 
bulk spin excitations, possible pinning of charge-density 
modulations (stripes), the local modulation of the su- 
perconducting order parameter, and large peaks in the 
density of states close to the Fermi level (as seen e.g. in 
STM). 

For smaller hole doping (i.e. in the insulating re- 
gime), the impurity effects on magnetism and trans- 
port have also been studied extensively. ^^'^"^'^^'^^ In 
La2-a;Sra;Cui_2Zn204 with x « 3%, co-doping with Zn 
leads to a re-appearance of ordered magnetism. The 
relatively high ordering temperatures suggest that this 
is not just the magnetism of impurity moments as de- 
scribed above for insulating magnets (which would be 
weak, as in Mg-doped TlCuCla^*). The continuous evo- 
lution of the ordering temperature with x and z also hints 
that the magnetism in the Zn and hole-doped material is 
smoothly connected to the commensurate antiferromag- 
netism of the parent compound, i.e., it is neither stripe- 
like nor strongly glassy. (Note that the freezing temper- 
ature in the spin-glass state of La2-xSra;Cu04 is much 
below the ordering temperature observed in Refs. 12,26. 
The suppression of glassy behavior upon Zn doping is 
also consistent with very recent /xSR and susceptibility 
data.^^) A natural conclusion is that the Zn impurities 
primarily "undo" the effect caused by the mobile holes on 
the commensurate antiferromagnetism. Resistivity mea- 
surements in La2-2:Sra;Cui_zZnz04 indicate a decreasing 
hole mobility with increasing Zn doping. Taken together, 
we conclude that Zn impurities tend to hinder the coher- 
ent motion of holes and thus prevent the magnetic order 
from being scrambled - this will be the central idea of our 
approach below. Other effects, namely the generation of 
effective staggered moments around Zn and vacancy di- 
lution of the magnetism (leading to percolation physics 
for z close to Zp), are clearly present, but are assumed to 
be subleading (at least for small impurity concentration). 

We note that an interesting dopant in La2Cu04 is 
Li (substituting for Cu) this provides both an extra 
hole and a vacancy in the CuO plane. Studies'^°'^^ of 
La2Cui_yLij^04 have shown that commensurate mag- 
netism survives up to y = 3 %. Transport data indicate 
that the holes introduced by Li doping remain localized 
also for larger y (in contrast to Sr doping) . Recent exper- 
imental studies show that charge and spin dynamics in 
La2Cui_j,Lij,04 appear to be glassy over a large range of 
doping. ■^^ Thus, the behavior of La2Cui_j,Lij^04 is differ- 
ent from the one of La2-:rSr,,.Cui_2Zn204 with x = z = y, 
the reason clearly being on the quantum chemistry side 
(e.g. the pinning potentials for holes being different for 
Li and Zn sites). Ref. 32 proposed an antiferromagnetic 
cluster state arising from Coulomb trapping of holes to 
explain some of the properties of La2Cui_j^Liy04, but 
more theoretical studies are clearly required. 
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C. Outline 

The remainder of the paper is organized as follows. In 
Sec. II wc will introduce the effective model describing 
the interaction between holes and spin waves. In Sec. Ill 
the hole and spin-wave Green's functions will be derived 
in the framework of the SCBA for the impurity-free sys- 
tem. Effects of the Zn impurities are discussed in Sec. IV, 
and Sec. V gives the expressions for magnetic properties 
within our approach. Finally, numerical results together 
with a comparison to experimental data are presented in 
Sec. VI. A discussion concludes the paper. 

We note that Korenblit et al.^^ have put forward a 
somewhat different explanation for the experiments of 
Hiicker et al.,^'^ based on the frustration model. ^'^ They 
argued that Zn impurities remove some of the frustrating 
bonds generated by mobile holes. However, their calcu- 
lations do not take into account the issue of reduced hole 
mobility due to impurity doping; we will further comment 
on their results towards the end of the paper. Impu- 
rity effects in cuprates have also been studied within the 
spin-fermion model. The authors do find a recovery of 
commensurate magnetism in a hole-doped situation upon 
adding impurities; however, the simplifications within the 
spin-fermion model (e.g. the absence of quantum spin 
waves) preclude a comparison with experiments. 



II. THE EFFECTIVE HAMILTONIAN 

In the following two sections we establish notation and 
- to keep the paper self-contained - summarize the mag- 
netic polaron model together with its SCBA treatment 
in the impurity-free situation, following Refs. 2,3,35. 

We start from the standard t-J model on a square lat- 
tice, H = Hf + Hj, and employ a representation using 
slave-fermion operators fi for spinless holes and Holstein- 
Primakoff bosons 6, for spin flips away from the Neel- 
ordered reference. This gives^'^^'^^ 

Ht = -^E/i 4(^+^1)-*' E /i/j+h-c, (1) 

(hi) {{hi}) 



(2) 



in standard notation. The first kinetic term in (1) de- 
scribes processes in which holes hop from one to the 
neighboring site creating or annihilating spin waves; the 
next-neighbor hopping term t' is allows for direct hole 
motion within one sublattice (an additional term of the 
form fifjbjbl is neglected). Values oi t/J = 3... 5, 
t' /t = —0.1 . . . — 0.2 are relevant for cuprate materials. In 
general, the Heisenberg term (2) contains an additional 
factor fifj fjfj which accounts for the loss of magnetic 
energy due to the hole doping. On the mean-field level 
we have fif- = 1 — ///i = 1 — <5, where S is the hole 
concentration. Consequently, the AF exchange coupling 
is renormalized according to J ^ (1 — 6)^ J, but for small 
hole doping we can neglect this effect. 



Fourier transforms of the b (/) are defined in the re- 
duced (magnetic) Brillouin zone, leading to b^, b^ {f^, 
f^) for the two sublattices A and B. The spin- wave part 
is diagonalizcd using a Bogoliubov transformation, intro- 
ducing new bosonic operators for spin waves, aq and /3q, 

with b^ = Uqaq + Wq/3^^q and 6^^ = v^a^ + itq/3^^, 
usual Bogoliubov parameters are given by 



The 



Ua = 



2Kc 



Vq = -Sgn(7q) 



2ko 



(3) 



with 7q = (cos^x + cosqy)/2 and Kq = (1 — 7q)^/^. Sub- 
stituting the new operators ctq, (3q into the Hamiltonian 
we arrive at the so-called spin-polaron model: 



q,k 

+y(-q, k - q)/^ V^q/3lq + h.c] 

k 

Hj = J2<«^^ + I3M- (4) 



Here, is the number of sites in each sublattice, and 
the coordination number is z = 4. All sums run over the 
magnetic Brillouin zone. The bare spin- wave energy is 



= ^/tq. The interaction vertex between holes and 



spin waves is given by F(q, k) = Mq7k-q + Wq7k and 
vanishes for q = or q = (tt, tt). As a consequence, not 
only the coupling between holes and long- wavelength spin 
fluctuations is important but also the coupling to short- 
wavelength spin fluctuations. 



III. GREEN'S FUNCTIONS AND BORN 
APPROXIMATION 

The hole Green's function in the absence of im- 
purities is defined by e^,(r,k) = -{T^f^{r)f^\0)), 
IJb,v = A, B, with the Fourier transform Q^^,{ivn,^ = 

/g''dr^^^(T, k)e*'^"'^, where Vn = {2n + l)'KT represents 
fermionic Matsubara frequencies. The spin-wave Green's 
function acquires the matrix form 



-(T,aq(r)«t(0)) -(T,«^(r)/?_q(0)) 
-{T,f3l^{T)alm -(r,/3l.q(T)/3_q(0)) 



The unperturbed spin-wave propagator is 

UJ,l -i^l^ 



Po^(iw„,q) 







(5) 



where cUn = 2iTnT, Introducing a self-energy matrix 
S(iw„,q), the solution of the Dyson equation reads 
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V{iuin,Cl) 



itj„ + w° + E22(«w„,q) -Ei2(iw„,q) 

-E2i(ia;„,q) -ia;„ + + Eii(iw„, q) 



(6) 



d(iw„,q) 
where 

d(ia;„,q) = [ia;„ + + E22(iWn, q)] - Wq - Eii(ia;„, q)] + Ei2(ia;„, q)E2i(ia;„, q). (7) 
The definition of the spin-wave propagators dictates that S22(«w„, q) = Eii(— ja;„, — q) and Si2(ia;„, q) = S2i(ia;„, q). 

A. Spin-wave self-energies 

The spin-wave self-energies, Eii(jcj„,q) and Ei2(ia;„,q), are calculated in the spirit of SCBA,''''' i.e., using the 
leading (bubble) diagram describing the decay of the spin fluctuations into a particle-hole pair with the full hole 
propagators, see e.g. Fig. 1 of Ref. 2. This approximation amounts to a summation of an infinite class of non-crossing 
diagrams; vertex corrections can be shown to be small for small doping^^ and will be neglected. Explicitly, the 
self-energies are given by 

Eii(ia;„,q) = (zt)^^ ^ ^ [T/(q, k)] ^ 5aa(«w„ k) ^bbC*!'™, k - q), (8) 

w„ k 

5^i2(iw„,q) = (2:t)^-^^-^^y(q,k)y(-q,k-q)^?AA(«w„-Fif„,k)^?BB(«t'n,k-q). (9) 

For simplicity and as in Ref. 3, in our numerics we will not perform a fully self-consistent calculation for the hole 
propagators, but instead use a trial form for the hole spectral density (derived from the SCBA solution of the single- 
hole problem) together with a rigid-band approximation. (A fully self-consistent SCBA calculation at finite hole 
concentration has been performed in Ref. 38 - it shows that at small hole concentrations the hole spectrum changes 
little with doping, i.e., the rigid-band picture is reliable.) 

I 



B. Hole spectrum 

The single- hole problem in the t-J model has been 
studied extensively using various techniques, includ- 
ing the SCBA within the spin-polaron model sketched 
above. '^^''^'''^^ It has been established that the spectral 
function for a single hole consists of a well-defined peak 
a low energies, /)coh(w,k), which corresponds to the co- 
herent motion of the dressed hole quasiparticle, and a 
broad incoherent part at higher energies, 

p{uj, k) = /9coh(t^, k) + Pincoh(w, k). (10) 

The coherent hole motion arises mainly from a com- 
bination of hopping and spin-flip processes and results 
in a narrow quasiparticle dispersion, e^! with a band- 
width of order 2 J and a dispersion minimum at momenta 
kmin = (if )if)- The quasiparticle weight, Zq, is re- 
duced from unity due to the dressing of the hole with 
spin fluctuation and scales with J/t in the physically rel- 
evant parameter regime. Therefore, for the coherent part 
we use^'^ 

Peoh(w,k) =Zo'5(a;-£k)e(2J-£k), (11) 

with the dispersion Ck = J(cos^ + cos^ ky). We have 
assumed that near the minima the dispersion is isotropic 



and can be approximated with £k « fc^/(2rneff), where 
k measures the distance from the kmin- Consequently, 
the effective mass of the hole scales as m~g ~ 2J. 
Here a remark is in order: In the t-J model without 
next-neighbor hopping, t' = Q, the dispersion near ki„i„ 
is known to be anisotropic, i.e., rather flat along the 
(0, 7r)-(7r, 0) direction.^^ However, photoemission results 
of Sr2Cu02Cl2 indicate an almost isotropic dispersion;^^ 
on the level of the t-J model a t' term is needed to capture 
this physics. '^^ Later on, it has been shown that a t-t'-J 
model reproduces salient photoemission features of differ- 
ent cuprate families.'*'* Therefore, we employ an isotropic 
dispersion in the following; we have also performed a few 
calculations with an anisotropic dispersion appropriate 
for t' = 0, leading to qualitatively similar results (the 6c 
in Fig. 1 below changes by less than 10%). 

The incoherent part of the spectrum arises from hole 
hopping inside the spin-polaron (or spin-bag) quasiparti- 
cle, with the characteristic energy scale t. A reasonable 
approximation is a momentum- independent constant: 

Pincoh(w) = e(w - 2 J) @{W + 2 J - to), (12) 

where W = 2zt is the incoherent bandwidth. 

Within a rigid-band approximation, the hole chemical 
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potential ^{6) is fixed by the hole density S through 

+ OC 

S=^Y1 / da;p(a; + /x,k)nF(a;), (13) 

-oo 

where p{lj, k) is the doping-independent hole spectral 
function (10), and n-p{uj) the Fermi function. For small 
hole doping {S <§; Zq) the chemical potential lies in 
the coherent band and at zero temperature has a value 
= ttSJ/Zq. 

Eqs. (6,8,9,13) completely describe the interaction be- 
tween hole quasiparticles and spin waves in the impurity- 
free host and have been used in Ref. 3 to investigate the 
destruction of magnetism in weakly doped cuprates. 



taken into account along the lines of Ref. 9. In Sec. VI B 
we will account for this dilution/percolation effect by 
simply multiplying our results with a magnetization re- 
duction factor taken from the studies of the diluted quan- 
tum Heisenberg antiferromagnet. 

Further, we will also neglect the physics of impurity- 
induced moments. We believe that magnetism aris- 
ing from these moments alone is too weak to explain the 
large observed Ncel temperatures, i.e., true bulk mag- 
netism is required. In a bulk-ordered phase, the effective 
impurity moments will order below Tn due to the mean 
field arising from the bulk order parameter this has 
been nicely seen^^ in La2-a;Sra;Cui_2Znj;04. However, 
the contribution of the impurity moments to the total 
order parameter is small. 



IV. THE EFFECT OF IMPURITIES 



B. Self-consistent T-matrix approach 



In this section we extend the theory to account for 

a finite density of non-magnetic impiirities in the sys- 
tem. As discussed in Sec. IB, we assiime that the vacan- 
cies primarily act as point-like potential scatterers for 
the mobile holes. Neglecting interference effects between 
spin-wave and impurity scattering processes for holes, we 
can capture the scattering off the impurities using a self- 
consistent T-matrix approach for the dressed hole quasi- 
particles, Eq. (10). 



A. Scattering potential 

Before proceeding a brief discussion of the quantum 
chemistry aspects is in order. Remarkably, there is no 
consensus on the description of Zn impurities within a 
one-band model for the CuO planes. Zn has a closed 
3d^^ configuration, suggesting that holes are expelled. 
This seems to be in agreement with cluster calculations 
in Ref. 46, which find the environment of a Zn atom to 
be similar to an undoped parent material. However, re- 
cent LDA calculations''^ have predicted a potential on 
the Zn site being repulsive for electrons (i.e. attractive 
for holes). In any case, the impurity potential is likely 
strong compared to the electron hopping. In what fol- 
lows we will model Zn as repulsive point-like scatterer for 
holes; for some parameter values we have also performed 
calculations with an attractive potential: the results will 
be mentioned below. (We note that Ref. 32 assumed a 
strongly attractive potential for holes associated to Li 
impurities in La2Cui_yLij,04 - if this is correct, and the 
potential of Zn is repulsive for holes, this would provide 
a reasonable explanation for the different properties of 
La2Cui_yLij^04 and La2-a;Sra;Cui_2Zn204.) 

In our calculations, we will neglect the effect of the im- 
purities on the spin waves. This piece of physics is well- 
studied:^' the magnetism is suppressed rather slowly 
with vacancy doping (compared to hole doping) . In prin- 
ciple, spin-wave scattering off the impurities could be 



For small impurity concentration, riimp <C 1, inter- 
ference effects between different impurities arc negligi- 
ble, and the self-consistent T-matrix approximation'*^ 
(SCTMA) is suitable to account for multiple scattering 
on a single impurity. Here we employ the SCTMA for 
spin-polaron quasiparticles, i.e., for holes which are al- 
ready dressed by the interaction with spin waves. Within 
SCTMA the full hole Green's function ^(iz/„,k) is given 
by 

g-\il^n,i^) = g-\iUn,-k) - S(iz/„,k), (14) 

where Q{ii'm'k.) is the hole Green's hmction correspond- 
ing to the spectral density (10), and the hole self-energy 
from impurity scattering is 



S(if„,k) = ni„ipTk,k(it'n)- 



(15) 



For point-like scatterers of strength Vq the T-matrix de- 
pends only on frequency: 



Vo 



1 - G{il^n) ' 



(16) 



and Q{ii'n) = N ^ G{ii^n, k) is the local hole Green's 
function. 

With the replacements ^ ^ 5 (14) in Eqs. (8), (9) 
and p ^ p = —'^mQ/iT in (13) we can now calculate 
how the impurities affect the spin-wave spectrum and 
consequently the magnetic properties of the AF host. In 
contrast to Ref. 2 we will not employ any further approx- 
imations. 



V. MAGNETIC PROPERTIES 

The spin-wave Green's functions defined in Sec. Ill al- 
low for a direct calculation of static magnetic properties. 
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A. Staggered magnetization 



The staggered magnetization in the Neel state de- 
creases with hole doping, 



m{S) — tuq — Am{S), 



(17) 



where the first term denotes the staggered moment in 2D 
AF reduced by ground-state spin-wave fluctuations only, 



1 



Too = 



1 

2N 



E 



1 



- 1 



(18) 



with K„ 



(1 



For a square lattice and spin- 



1/2 we have the well-known spin- wave result too ~ 0.303 
(which is in excellent agreement with toq = 0.3070(3) 
from quantum Monte Carlo simulations^"). The effect 
of doped holes on the AF order is included through the 
second term in (17), given by:"^^ 



q 



7q(aq/3- 



(19) 



The bosonic occupation numbers in (19), being zero in 
the undoped case at T = 0, can be calculated from 
the spin- wave Green's functions as {a^^a^) — ~'Dii{t — 
0~,q), (aq/3_q) = —'D2i{T = 0~,q). The phase boundary 
of the antiferromagnetic phase is given by to ((5c) = 0. 

Notice that m{S) (17) refers to the magnetization per 
spin (within the spin- wave approximation). To obtain 
the magnetization per site, as measured by bulk probes, 
m{S) has to be multiplied with the factor [l — d — nimp); 
this does not change the location of the phase bound- 
ary, dc- Of course, this mean-field factor only crudely 
accounts for the dilution of the magnetic lattice; in par- 
ticular, percolation physics cannot be captured. We will 
comment on this in Sec. VI B 2 below. 
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FIG. 1: (color online) Doping dependence of the staggered 
magnetization for "mobile" (i.e. dispersive) holes (in the ab- 
sence of impurities) for different values of J at T = 0. (Typ- 
ical parameters for cuprates obey t/ J — 3 . . .5.) The weight 
of the dispersive quasiparticle peak in the hole spectrum is 
Zo = 0.5 J/t; the lines are guides to the eye only. The data 
are very similar to the ones of Ref. 3. 



VI. NUMERICAL RESULTS 

In this section we will present our numerical results 
for the doping dependence of the AF order parameter 
with and without impurities. Note that all energies are 
given in units of t. Most calculations are done for a 
two-dimensional system at zero temperature;^^ the finite- 
temperature results are obtained using a finite interlayer 
coupling. 



B. Neel temperature 

The above calculations, performed at finite tempera- 
ture, allow to determine the Neel temperature, Tn, from 
m{S, T = Tn) = 0. This requires a weak AF coupling, J±, 
along the direction perpendicular to the CuO plane. In 
high-Tc cuprates this inter-layer AF coupling is typically 
J_L ~ (10~^ — 10~^) J. As a consequence, the spin-wave 
dispersion has a quasi-2D form which is given by 

(20) 

where 

7q - 7q + e7q • (21) 

Here 7^ = cos Qz , and the ratio between inter- and intra- 
layer coupling is denoted by e = J±/{2J). 

A drawback of the spin-wave approach is that no sen- 
sible critical behavior is obtained, i.e., the magnetization 
vanishes linearly as T — s- Tn (the order parameter critical 
exponent is P = 1). 



A. Hole-doped antiferromagnet without impurities 

First we will discuss the case of the doped AF host 
without impurities and compare our results with pub- 
lished ones.^'"^ Extending these calculations, we will con- 
sider different forms of the hole spectrum to study the 
influence of both hole dispersion and hole coherence on 
the host magnetism. 



1. "Mobile" holes 

The situation of "mobile" holes refers to a hole spec- 
trum function as obtained from solving the single-hole 
problem in the t-J model, i.e., the hole spectral func- 
tion consists of both a coherent, dispersive (11) and 
an incoherent, localized (12) part. In Fig. 1 the zero- 
temperature doping dependence of the staggered magne- 
tization is shown for different values of the coupling J 
and fixed quasiparticle weight Zg (taken from a single- 
impurity calculation'^^'^^''"'). 
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FIG. 2: (color online) Doping dependence of the staggered 
magnetization now for "localized" holes, i.e., in the absence of 
a coherent quasiparticle peak in the hole spectrum {Zq — 0). 
The calculations have been done at temperature T — 0.005, 
for fixed J± = 0.01 J. 



For J — 0.20 we obtain the similar critical hole concen- 
tration as Belkasri and Richard,'^ Sc ~ 2.7%. Our result 
for J = 0.25 is a bit smaller than one derived by Khal- 
iuUin and Horsch,^ (5c ~ 4%, due to the fact that (as a 
further approximation) they have taken into account only 
the spin- wave self-energies with small momenta. The dis- 
appearance of the AF order at hole concentrations Sc of 
a few percent is in good agreement with experiments on 
both La2-a;Sra;Cu04 and YBa2Cu306-i-5. 



2. "Localized" (incoherent) holes 

Now we turn to a limiting case where the holes are lo- 
calized in a certain extended region and only move inco- 
herently (e.g. around a pinning center). With the above 
form of the spectral function (10) this corresponds to the 
quasiparticle weight being zero, Zq = 0. Thus, we use 
p(aj,k) = pincoh(w) (12). We study this toy situation to 
illustrate the strong influence of hole localization on the 
magnetic properties. 

One expects the incoherent motion of holes is less de- 
structive for the AF host in comparison with the case 
of "mobile" holes. In Fig. 2 we show the doping de- 
pendence of the staggered magnetization for "localized" 
holes and different values of AF couphng J. {t/J is still 
a relevant parameter here because the bandwidth of the 
incoherent motion is controlled by t. Also note that re- 
sults are shown for finite temperatures and finite cou- 
pling J± implying, e.g., that the value of the staggered 
magnetization m(0) differs from the T = value mo. 
However, this difference is rather small since the effects 
of T > and J± > tend to compensate.) The crit- 
ical hole concentrations 5c here are significantly larger 
than for dispersive holes, but still much smaller than for 
strictly static vacancies (where Sc equals the percolation 



threshold Zp = 40.5 %). The reason is of course that even 
our spatially localized holes scramble the magnetic back- 
ground in their vicinity. We call the reader's attention 
to the analysis of hole mobilities in Ref. 12, which con- 
cluded that "completely localized" holes (i.e. with zero 
mobility) still destroy the host magnetism three times 
faster than static vacancies - this gives 5c ~ 13 % in rea- 
sonable agreement with our data in Fig. 2. 

In Fig. 3 the doping dependence of the Neel tempera- 
ture, Tn, for the "localized" holes is depicted for the same 
set of parameters as in Fig. 2. As also known for the case 
of doping with mobile holes, ^'^ the initial suppression of 
the Neel temperature with doping is larger than the one 
of the staggered magnetization. 

Comparing the results for "mobile" and "localized" 
holes, Figs. 1 and 2, one can easily conclude that inco- 
herent motion of holes is less destructive than coherent 
one. Therefore, any mechanism which modifies the hole 
spectral function in such a way that suppresses coherent 
and enhances incoherent motion will lead to a reinforce- 
ment of commensurate antiferromagnetism. (Note that 
the limit t — > 0, corresponding to immobile holes and 
percolative destruction of magnetism, is not described 
within the current approximation, but can be captured 
as in Sec. VI B 2 below.) 



B. Hole-doped antiferromagnet with non- magnetic 
impurities 

Let us finally turn to our main novel results, namely 
the magnetic properties of the system with simultaneous 
hole and impurity doping. Here, the spin- wave properties 
are calculated using a hole spectrum function correspond- 
ing to "mobile" holes, i.e., consisting of a quasiparticle 
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FIG. 3: (color online) Doping dependence of the Neel temper- 
ature for "localized" holes {Zo = 0), with fixed J± = 0.01 J. 
(Typical Tn values for undoped cuprates are a few hundred 
kelvins, e.g., La2Cu04 has Tn « 300 K. The hopping in- 
tegral t in the t-J model for cuprates is 300 . . . 400 meV 
^ 3000 . . . 4000 K.) 
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peak and an incoherent background, but subjected to im- 
purity scattering (which re-distributes weight in the hole 
spectrum) . 

Typical results for the AF order parameter as function 
of the hole doping are shown in Fig. 4. Clearly, Zn doping 
increases the critical hole concentration Sc- The differ- 
ence between the Sc in a case with and without impurities 
is of order of one percent which is in reasonable agree- 
ment with experiments^^ - this indicates that we indeed 
capture the main effect of Zn doping in La2-xSra;Cu04. 
Notice that 5c increases only weakly as function of the 
scattering potential Vq for Vb > 4. We have performed 
calculations for different values of J, with qualitatively 
similar results (all values increase for larger J, see 
Fig. 1). In addition, we found that small to moderate 
negative Vq has qualitatively the same effect as positive 
Vq (strong negative Vq leads to hole bound states, where 
holes get trapped, and consequently a larger increase of 
6c through impurity doping). 

1. Less coherent holes 



-.- Vo 


= 


— Vo 


= 1.0 


* Vo 


= 2.0 


-T- Vo 


= 10.0 
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S 

FIG. 4: (color online) Doping dependence of the stag- 
gered magnetization for simultaneous hole and impurity dop- 
ing. The different curves correspond to different impurity 
strengths Vo, the impurity concentration is nimp = 15%, the 
Vo ~ curve is the impurity-free case. The other parameters 
are J = 0.20 and Zo = 0.5 J/t, T = 0. 



Within our calculations we can clearly identify the rea- 
son for the antiferromagnetism being more robust in the 
presence of impurities. In the self-consistent T-matrix 
approach of Sec. IV the main effect of the impurity scat- 
tering is to transfer spectral weight, i.e., to reduce the 
weight of the coherent dispersive quasiparticle peak. We 
can quantify that: for ni,„p = 15 % and Vq = 2 the quasi- 
particle weight Zq is reduced by roughly 10%. A toy 
calculation with this reduced weight (but in the absence 
of impurities) shows that this loss of quasiparticle weight 
accounts for 3/4 of the change in m between the clean 
and impurity-doped cases. (The remaining 1/4 arises 
from a change in the quasiparticle dispersion and from 
re-shuffling spectral weight at higher energies.) 

Thus, the loss of coherence of the hole quasiparticle due 
to impurity scattering is the main source of the reinforce- 
ment of magnetism. We can interpret this in real space: 
The coherent motion of the hole quasiparticle in the an- 
tiferromagnetic background corresponds to the emission 
of a spin wave in each hopping step, with a strong sup- 
pression of the order parameter. In contrast, during the 
incoherent motion within the quasiparticle the hole also 
absorbs spin waves (due to self-retracing paths), and the 
suppression of antiferromagnetism is less severe. 

2. Spin-wave scattering and percolation 

Our approach has so far neglected the dilution effect 
of the vacancies on the magnetic, i.e., spin- wave, part of 
the model. On the mean-field level, dilution results in 
a multiplicative factor of (1 — ^ — ^imp) for m to obtain 
the staggered magnetization per site. However, impu- 
rity scattering of spin waves is still not included, and 
our approach falls short of capturing percolation physics. 



While the impurity effect on spin waves could be treated 
in principle,^ we choose a different route here: As dilu- 
tion of a well-ordered antiferromagnet has been studied 
extensively,^" we can simply multiply our magnetization 
results with a renormalization factor P(nimp) which rep- 
resents the (normalized) magnetization of a 2D square 
lattice quantum Heisenberg magnet diluted with Tiimp 
vacancies - this quantity can be taken from quantum 
Monte Carlo simulations. Fig. 18 of Ref. 10. With this 
we reproduce a vanishing magnetization at the percola- 
tion threshold.^'* 



3. Re-appearance of magnetism 

To emphasize the fact that impurities help the AF or- 
der to recover we show in Fig. 5 the staggered magnetiza- 
tion m at zero temperature as a function of hole doping 
with and without impurities - here m is not normalized 
to its value at 5 = 0. Another way to show the recovery 
of long-range order in the presence of impurities is pre- 
sented in Fig. 6, where order parameter is plotted as a 
function of impurity concentration for two different hole 
dopings. For 5 = 2.7% the AF order sets in at finite 
impurity concentration, riimp « 2.5 %. The presence of 
impurities causes the re-appearance of the AF order for 
hole concentrations greater than S w 2.6%, and the or- 
der parameter depends non-monotonically on the hole 
doping - both results are in good agreement with the ex- 
periments done on La2-a;Sr2;Cui_^Zn204 by Hiicker et 
al}^ 

For a detailed comparison with experiments a few re- 
marks are in order: (i) The present calculation is strictly 
valid only for small impurity concentration, riimp <ti 1, 
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FIG. 5: (color online) Re-appearance of the AF ordering 
through impurities: The staggered magnetization, corrected 
by the percolation factor P(ninip) (see text), as function of 
hole doping for zero and 15 % impurity concentration. The 
parameters are J = 0.2, quasiparticle weight Zo = 0.5 J/t, 
impurity potential Vb = 2.0, T = 0. 



when the interaction between impurities is neghgible. (ii) 
Our simpHfied treatment of magnetic dilution, utihzing 
numerical results for the diluted insulating Heisenberg 
antiferromagnet, only partially captures the interplay of 
dilution and quantum fluctuations. In particular, in the 
presence of vacancies and a small amount of holes mag- 
netic order does not need to persist up to the percola- 
tion threshold.^"* Thus, our data at finite hole doping 
become unreliable for riimp > 25 %. (iii) Corrections be- 
yond spin-wave theory are required for the description 
of critical behavior in the vicinity of Sc', we expect m 
to vanish with an infinite slope (i.e. the order parameter 
exponent (3 being smaller than unity), (iv) Experiments 
usually measure the ordering temperature Tn, not the 
zero-temperature order parameter m. Although the two 
behave in a qualitatively similar fashion, Tn decreases 
initially faster with doping than m. Also, for a complete 
picture it may be necessary to include the interlayer cou- 
pling beyond mean- field. The detailed behavior of Tn 
with both hole and impurity doping, taking into account 
percolation-type physics, is beyond the scope of this pa- 
per. 

In Fig. 7 we show the evolution of the critical hole 
doping level with the impurity concentration. Inter- 
estingly (and in agreement with experiment), the varia- 
tion of Sc is rather small, i.e., vacancies cannot be used 
to shift Sc to arbitrarily large values. 

Last not least we briefly comment on the calculation 
by Korenblit et al.,^^ which provides a different explana- 
tion for the experiments of Ref. 12. They present scal- 
ing arguments based on the frustration model, with the 
key idea that Zn removes ferromagnetic bonds created 
by hole doping. This simple theory predicts that Sc can 
be increased to large values by introducing Zn impurities 



FIG. 6: (color online) Same as Fig. 5, but now showing m as 
function of the impurity concentration for different fixed hole 
dopings. m vanishes at the percolation threshold, nimp = 
40.5%. 



this is in contradiction to the experiment. Thus, their 



y cannot apply to 



simple scaling form S — > S{1 — riimp 
larger S and Tiimp- Also, the paper makes no reference 
to the motion and the mobility of holes - we think that 
this aspect is crucial for the description of hole-doped 
cuprates. 



VII. CONCLUSIONS 

Motivated by experiments on Zn-doped 
La2_a;Sra;Cu04, we have discussed the interplay of 
mobile holes and static non-magnetic impurities in 
lightly doped antiferromagnetic Mott insulators. Based 
on a self-consistent spin-wave and T-matrix calculations 



<5c 




FIG. 7: (color online) Critical hole doping level, Sc, as func- 
tion of the impurity concentration nimp for different values of 
the scattering potential. The parameters are J = 0.2, quasi- 
particle weight Zo = 0.5 J/t, T = 0. 
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wc have shown that commensurate long-range order 
is more robust against hole doping in the presence of 
non-magnetic impurities as compared to the clean case. 

We believe that the most naive picture, namely that 
Zn impurities trap holes, is not appropriate. This can 
already be seen from the experimental data: at a hole 
doping of 2.3 % the Neel temperature depends strongly 
and non-monotonically on the impurity concentration for 
4 % < riiiiip < 15 %. Instead, wc propose that impurities 
mainly reduce the hole mobility and suppress the coher- 
ent part of the hole motion (in favor of an incoherent 
background) . In turn, this suppresses the spin- wave soft- 
ening caused by coherent hole motion. In the case of 
Ni doping^'^ this effect is corroborated by the ordering 
tendencies of the spin-1 impurity moments. 
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